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An inverse Compton scattering (ICS) model of pulsar emission: 

III. polarization 

R. X. J. F. Liu2'3, J. L. Han^'^, G. J. Qiao^'^'S 

ABSTRACT 

Qiao and his collaborators recently proposed an inverse Compton Scattering 
(ICS) model to explain radio emission of pulsars. In this paper, we investigate the 
polarization properties of pulsar emission in the model. First of all, using the lower 
frequency approximation, we derived the analytical amplitude of inverse Compton 
scattered wave of a single electron in strong magnetic field. We found that the 
out-going radio emission of a single relativistic electron scattering off the "low 
frequency waves" produced by gap-sparking should be linearly polarized and have no 
circular polarization at all. However, considering the coherency of the emission from 
a bunch of electrons, we found that the out-going radiation from the inner part of 
emission beam, i.e., that from the lower emission altitudes, prefers to have circular 
polarization. Computer simulations show that the polarization properties, such as the 
sense reversal of circular polarization near the pulse center, S-shape of position angle 
swing of the linear polarization, strong linear polarization in conal components, can be 
reproduced in the ICS model. 

Subject Headings: polarization — Pulsars: general — Radiation mechanisms: 
non-thermal 



1. Introduction 

The outstanding polarization properties of radio pulsars are keys to understand the 
magnetospheric structures of the neutron stars and the unknown emission mechanisms. Not only 
linear but also circular polarization has been detected from most pulsars. Soon after the discovery 
of pulsars, Radhakrishnan & Cook (1969) proposed the rotation vector model to interpret the 
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rotating position angle of linear polarization, which has been widely accepted for practical reasons, 
regardless of the emission mechanism. Obviously, the linear polarization is probably related to the 
structure of the strong dipole field above magnetic poles. Large amount of polarization data have 
been accumulated (e.g. Gould & Lyne 1998; Manchester, Han &; Qiao 1998; Rankin, Stinebring 
& Weisberg 1989; Weisberg et al. 1999). However, many observed polarization features can not 
be well explained, for example, the orthogonal position angles of linear polarization observed from 
individual pulses (e.g. Stinebring et al. 1984a,b), diverse circular polarization (e.g. Han et al. 
1998), the different polarization characteristics of core and conal emission (Rankin 1993; Lyne 
Sz Manchester 1988). Though there were many theoretical efforts (e.g. Gil & Snakowski 1990), 
no consensus has ever reached on how pulsar polarization was generated (Radhakrishnan 1992; 
Melrose 1995). For example, the observed circular polarization in pulsar radio emission could 
either be converted from the linear polarization in the pulsar magnetosphere via the propagation 
effect (Melrose 1979, Cheng & Rudcrman 1979; von Hoensbroech & Lesch 1999) or plasma process 
(e.g. Kazbegi et al. 1991), or origin from emission process intrinsically (Radhakrishnan &; Rankin 
1990; Michel 1987; Gangadhara 1997). 

Recently, Qiao and his collaborators proposed an inverse Compton scattering (ICS) model of 
pulsar emission (Qiao 1988; Qiao & Lin 1998; Qiao et al. 1999). In the model, radio emission is 
the result of the secondary particles scattering off "low-frequency waves" . The waves are assumed 
to be produced by either the breaking down of the vacuum polar gap, or other sorts of short time 
oscillations or micro-instabilities (Bjornsson 1996) . The gap-sparking may exist above polar cap, 
as shown by new observations by Deshpande & Rankin (1999) and by Vivekanand & Joshi (1999). 
If so, there must be low-frequency waves emitted by such sparking. The time scale of a sparking at 
a certain point is about 10~^ s or even shorter, the out-flowing plasma produced by the cascades 
of the sparking should form miniflux tubes with a radius of ~ 10 m and a length of ~ 100 m, 
and hence be highly inhomogenous in space. The plasma density between the tubes should be 
sufficiently low for radio wave to propagate, as if in a vacuum. Therefore, we assume that the 
pulsar magnetosphere is transparent for radio emission. That is to say, electromagnetic waves of 
observed emission can pass through the magnetosphere freely, and the so-called "low-frequency 
wave" can propagate through the emission region above the polar gap (sec Qiao & Lin 1998). 
The success of the ICS model (Qiao & Lin 1998, Qiao et al. 1999) is that the core and the conal 
emission beams can be naturally explained in this model, while the core component comes from 
the nearest region to neutron-star surface, the inner cone is farther, and the out cone is farthest. 
Because of different heights of these emission components, the retardation and aberration effects 
cause the components apparently shift spatially, and could produce the position angle jumps in 
the integrated pulse profiles in some cases (e.g. Xu et al. 1997). 

In this paper, we investigate the polarization properties of the ICS model for radio pulsars. 
First of all, we presented the polarization characteristics of scattered radio emission from a 
single particle in Sect. 2, half from theoretical work, half from numerical calculations. In the 
low frequency limit, deduction of the amplitude of outgoing radio waves from the ICS process 
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in classical electrodynamics (CED) was presented in Appendix. In Sect. 3, we showed that the 

circular polarization can be the result of coherency of the outgoing waves. Considering the 
geometry above the polar cap regions, we numerically simulated the ICS process, and found that 
the circular polarization with possible sense reversal is preferably found from the emission in the 
beam center, while stronger linear polarization preferably from the outer part of the beam. Some 
related issues will be discussed in Sect. 4. 



2. Polarization features of the out-going scattered waves from a single electron 

The polarization features of scattered emission by relativistic electrons in the strong magnetic 
field has not been discussed in literature. In this section we are concerned with the general 
polarization properties of the scattered photons by a single relativistic electron onto the lower 
frequency waves. Previously, the quantum electrodynamics (QED) results of Compton scattering 
cross sections were first derived by Herold (1979) for electrons with ground initial and final 
states. Then, the QED cross sections for various initial and final electron states were calculated by 
Melrose & Parle (1983) and Daugherty &; Harding (1986). Finally, Bussard, Alexander &; Mesaros 
(1986) presented the QED results for electrons with arbitrary initial and final states. Xia et al. 
(1985) calculated for the first time the differential cross section of inverse Compton scattering, 
but no polarization features were discussed. 

Based on the results of cross sections of various polarization vectors derived by Herold (1979), 
we used the Lorentz transformation in the magnetic field direction to calculate the differential 
cross section of inverse Compton scattering of orthogonal polarization pairs, an approach similar 
to that in Xia et al. (1985). Hereafter we denote the angular frequency of an incident photon as 
Win, the strength of magnetic field as B, and its cyclotron frequency as = and the angular 
frequency of the scattered wave or photon as Wout- The geometry of the ICS process is shown in 
Fig. 1. The relativistic electron is moving along the field line (on the z axis) and has a Lorentz 
factor 7 (and = 1 — 1/7'^ in the equations) which we suppose to be about ~ 10^ to ^ 10^ in 
the pulsar magnetosphere. From the kinematic restrictions for the scattering process in strong 
magnetic fields, i.e. the momentum conservation in the magnetic field direction and the energy 
conservation of the process, the frequency of the scattered photons is (note that "sin0" in Eq.(15c) 



where / = 7(1 — P cos 9out)n^c'^ + (1 ~ cos cos ^out)^"^!!! and g = ^Win sin^ 0out[27(l — 

/3cos ^in)mc^ + sin^ ^in^f^in]//^- When <S rrieC^, all above three equations can be simplified as 



of Xia et al. 1985 should read as "sin' 





(1) 
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being 

1 - /? cos 6'in 
1 - /3 COS 6'out 

According to the definition given by Saikia (1988), we derived the Stokes parameters (similar 
to Eq.(6) below) of scattered wave with the geometry plotted in Fig. 1. We numerically calculated 
the polarization properties, using Herold's results, for radio frequency regime and high energy 
bands. We found that the polarization properties of the scattered photons are different. 

(1). When Win "^c and both and Wout are in radio band, the scattered photons are 
completely linearly polarized, and its polarization position angle (PA) is definitely in the co-plane 
of the out-going photon direction and the magnetic field (see geometry in Fig. 1), in spite of the 
polarization state of incident photons with one exception for the very specific case in which the 
incident photon is absolutely linear polarized and its polarization vector is exactly perpendicular 
to the magnetic field line. 

We derived the amplitude of the scattered waves from a single electron with the classical 
electrodynamics (CED), which we will use in next sections for emission coherence. The electron 
can be considered as a particle since the magnetic field B <^ Bg, here Bg is the critical magnetic 
field Bg = {vn?c^)/{e?i) = 4.413 x 10^^ G, and the incident and the scattered photons can be 
considered as waves since the frequency of the incident wave in the electron-rest-frame is much 
smaller than the resonant frequency, i.e., Win7(l — f3cos9\^) ^ Uc- We derived (see Appendix) that 
the electric field Es(t) of the scattered wave is 

p. (f\ - ''^ p cos 7] sin gjn sin Op^t ^i(k.D-a,outt) ^out /oX 

where k the wave-vector, D the position vector from scattering electron to observer, Tc the electron 
classical radius. The incident wave has electric field amplitude Eq, and its electric vector has an 
angle r/ with respect to e™. The wave is scattered to be the out-going wave with only completely 
linear polarization, the polarization vector of which (ej"^') is in the plane of out-gong wave and 
magnetic field, in consistent with the results from the QED calculation shown above. We note 
that Chou & Chen (1990) have calculated Stokes parameters of Thomson scattered wave in 
strong magnetic fields and found that the scattered radiation is always linearly polarized for any 
polarized incident wave. Their results serve as an independent check of the specific case of 7 = 1 
above. 

In 3-D, emission from a single electron is going out in a micro-cone around the magnetic line, 
the cross-section profile is shown in Fig. 2, with the maximum at about ^?out = — O.6/7 (i.e. 
0.3° when 7 = 100. See Fig. 2), as deduced from Eq. (A9). At this angular radius, the frequency 
of the scattered waves should be (from Eq. 2) 

Wout = 1-5 a;in7^(l - l5cosBiri)- (2a) 



When a line-of-sight goes across the micro-cone, one can detect an "S" shape variation of 
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polarization angle. Note that this cross section of the micro-cone does not have a gaussian-shape. 
All those polarization properties do not vary with frequency until u-m is about 10^^ s~^. 

(2) . For resonant scattering at high energy bands (for B = lO^^G, ujin = 1.3 10^^ s~^), the 
scattered emission (wout = 1-7 10^^ s^^) has a filled gaussian-like shape with the total emission 
peak in the direction of magnetic field (sec Fig. 3), rather than a micro-cone. Circular polarization 
is 100% near the center, but decreases to zero when the angular radius ^out is about I/7, and it 
changes the sense and gets almost 100% again when ^out = 1-5/7. 1^ contrast, there is no linear 
polarization near the center and the edge, but it peaks up to almost 100% when ^out = 1/7- The 
position angle of linear polarization of the scattered photon is perpendicular to the co-plane of 
out-going photon and the magnetic field, different from the case in radio band. 

(3) . When win is about the resonant frequencies, i.e. higher than about 10-*^^ or lower 
than 10^*^ s^^ as seen in Fig. 4, the total intensity of the scattered photons increases towards the 
resonant frequency. However, this increased total intensity is mainly contributed from circularly 
polarized emission, which is shown in the upper plots of Fig. 4 from the fraction variations of 
circular polarization with frequency at given directions of out-going photons. 

What we considered in the second and third cases is at high energy regime with parameters 
B = 10^^ G and 7 = 100 in this section. However, polarization properties of scattered emission 
we presented above should be the same for the situation in the resonant region in the outer 
magnetosphere of pulsars which we will not considered in this paper. 

Note that all these conclusions for an electron are the same for a positron because the 
amplitude of the scattered waves rests on the incident waves, having nothing to do with the charge 
sign of a particle. 

3. The coherent ICS process in pulsar magnetosphere 

Strong polarization is the most outstanding feature of pulsar emission. Circular polarization 
has been detected from pulsars, and sometimes has sense-reversals across pulse profile (Han et al. 
1998). In this section, we consider the pulse polarization in the frame of the ICS model. 

As we see from the last section, scattered emission in radio band from a single electron is 
completely linearly polarized, and the position angle follows an "S" shape for a sweep of the line 
of sight. The scattered waves from electrons in a close group are coherent, which can produce the 
circular polarization detected from pulsars. The amplitude of coherent emission can be calculated 
by using Eq. (3) for radio emission, so that we can avoid the much more complicated procedures 
by using coherent states in quantum electrodynamics. For example, if two linearly polarized plane 
waves El and E2 with same wave vector have a phase difference 5 and an angle of k between their 
polarization vectors, the circular polarization percentage of the coherently superimposed emission 
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is 

V 2E1E2 sin K sin 5 

I ~ El + El + 2E1E2 cos K cos 5 

For El = E2, Eq.(4a) becomes 

V/I = sinKsmS/{l + cosKCOsS). (46) 

This gives considerably circular polarization, as long as the values of | sin«;| and \smS\ are not too 
small, li K = S = 7r/2, it will be totally circularly polarized. 

In the ICS model of radio emission of pulsars, the incident waves were generated by a 
sparking due to the breaking down of inner gap. The waves travel through a vacuum-like pulsar 

magnetosphere with a small filling factor of pair plasma (i.e., inhomogenous spatially). They 
encounter a bunch of relativistic particles that come from another cascade (in the other side) above 
the polar cap, and arc upscattcrcd coherently (Fig. 5) to produce the observed radio emission. In 
a given observational direction, coherently superpose scattered waves from a bunch of particles, 
which is seen from the transient "mini-beam" below, can be almost completely linearly polarized 
or some times have very strong circular polarization dependent on geometry. This is mainly the 
result of the coherency and the non-random distributions of the 5 and k in Eq.(4). However, the 
observed integrated pulse profiles of pulsars, is the incoherent sum of many samplings of these 
mini-beams. In the field of view of a given line of sight, if the sparkings distribute symetrically, the 
circular polarizations will cancel each other, and no circular polarization left finally. Nevertheless, 
there is remarkable circular polarization in mean pulses if the sparking distribution is asymmetric. 

In this section, we first describe the coherent ICS process of a bunch of particles, then we 
present computer simulations of their transient beam. Finally we simulate the mean pulse profiles, 
which have linear and circular polarization, S-shaped position angle swing and unpolarized 
emission as well. 

3.1. The scattered waves from a bunch of particles 

The scattered emission from a single particles is in a micro-cone of about O.6/7 (see Fig. 2 
and Sect. 2), hence only emission from a small area is visible to a given line-of-sight hq. Since a 
particle is moving in a direction along the magnetic field hb, the emission from these particles 
satisfying • no ~ cos(0.6/7) can be received by an observer. Using Eq.(2a) (see also Eq.(6) 
of Qiao & Lin 1998), combined with Eqs.(ll) and (13) of Qiao Sz Lin (1998) which describe the 
scattering geometry and how the Lorentz factor of the particles changes along the field lines, one 
can find three emission heights for a given frequency. At a certain observing frequency, an observer 
cannot receive emission at the same time from all zones at three heights produced by a bunch of 



(4a) 
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particles, but just see a part of one zone 0. 

The upper-scattered waves from particles in such a small area can be accumulated coherently 
for some good reasons. First of all, the low frequency incident waves that were generated by one 
sparking should have a phase coherence when they encounter the bunch of particles, even if with 
a small phase dispersion due to travelling. Secondly, the scattered waves at one given frequency 
from these particles in a small visible ring-like area would have a good phase coherence as wel]|^. 
The complex electric field amplitude contributed from each particle[see Eq.(3)] reads 

COSr/sine^in r-/'^inj5 ^out „ , , x-, out (r, \ 

Es = C exp z( — R R • no + 0o)] , (3a) 

R c c 

here C is a constant for a given 7. Vectors R = SA and no are illustrated in Fig. 5, and (pQ is 
the initial phase of the incident waves when it was generated by a sparking. Note that (pQ might 
be random for different sparkings. When the incident waves come from one sparking, we can 
integrate the complex amplitudes (rather than the power) from Eq. (3a) of scattered waves from 
all visible particles. The complex electric field amplitude of the total scattered wave is, 



E = / B^riedV. (5) 

J nB-no~cos(0.6/7) 

Here, rig is the number density of particles in the bunch. The Stokes parameters of this 
accumulated emission are (Saikia 1988) 



I = i[(eio-E)*(eio-E) + (e2o-E)*(e2o-E)], 

Q = i[(eio-E)*(eio-E)-(e2o-E)*(e2o-E)], 

U = Re[(eio-E)*(e2o-E)], 

V = Im[(eio-E)*(e2o-E)]. 

Two orthogonal polarization vectors above were defined to be e2o = no x S7/|no x ft\, and 
Gio = ^20 X no, here $7 is the direction of rotation axis (Fig. 5). 



(6) 



3.2. Transient beam from a bunch of particles: simulation 

We now simulate the coherent ICS process of a bunch of particles, and will present its 2-D 
emission feature snap-shot. 

The particles produced by the sparkings are moving out along the extremely strong magnetic 
fields. The radius of a sparking spot was assumed to have the same scale as the inner gap height 



^ Actually, one just see a ring-like area about n_B ■ no ~ 003(0.6/7), where (1) the scattered emission is at a given 
observational frequency; (2) the power reaches its maximum. 

® We understand that if this area is much less than the wavelength of incident waves, there would be coherence. 
Otherwise, for example, "mini-beam" from the outer cone region we discuss below is much larger, and can't be treated 
coherently. 
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(Gil 1998) , which Deshpande & Rankin (1999) estimated to be about 10 meters. We took this 
value in the following simulations, but we have found that changing of this dimensional size does 
not affect our results. We also assumed that the number density of particles rig in a bunch has 
the maximum in the center and decline towards the edge. We took a gaussian distribution for 
convenience. This number density will be the natural weight upon the contribution in Eq.(5) from 
different visible part of a ring-like area at a given frequency. The wave-phase of this contribution 
depends on the location of this visible part relative to a sparking spot. 

The ICS process by one bunch of particles moving along given field lines forms three 
"mini-beams" of a few degrees in width. They were produced at emission heights of core and two 
cones, respectively. We assumed that the particles in the bunch have the same Lorentz factor, and 
found that the polarization features are quite different between these mini-beams(see Fig. 6). 

As seen in Fig. 6, circular polarization is very strong, even up to 100%, in the core minibeam, 
as shown by the open and filled circles. Emission from the minibeam of inner cone is much less 
circularly polarized. However, we could not calculate the beam feature of the outer cone using 
coherent superposition method as in the case of core and inner cone, since the particles at this 
larger emission height have much smaller Lorentz factors, and therefore the visible ring-like area 
(see Fig. 5) is so large that the upscattered emission is no longer coherent. 

We also made simulations for a bunch of particles with different Lorentz factors, but the 
polarization features shown in Fig. 6 do not change significantly. This does not surprise us since 
the global polarization pattern is mainly determined by field geometry. 

When the line of sight sweeps across a minibeam, we can see a transient "subpulse", as 
shown in Fig. 7. Note that the width of "transient" subpulses is about 0.6°, much larger than 
1/7 ~ 0.04°. When the line of sight sweeps across the center of a core or inner conal minibeam, 
the circular polarization will experience a central sense reversal, or else it will be dominated by 
one sense, either left hand or right hand according to its traversal relative to the minibeam. 

The position angles at a given longitude of transient "subpulses" have diverse values around 
the projection of the magnetic field (see Fig. 6 and Fig. 7). The variation range of position angles 
is larger for core emission, but smaller for conal beam. When many such "subpulses" from one 
minibeam is summed up, the mean position angle at the given longitude will be averaged to be 
the central value, which is determined by the projection of magnetic field lines. Therefore in 
our model, the mean position angle of pulsar emission beam naturally reflects the geometry of 
magnetic field lines. As we will show in our simulation the position angle of mean pulse profiles 
should have "S" shapes (see section 3.3). 

Note that these transient "subpulses" are not the subpulse we observe from pulsars. The 
time scale for the transient beam|] to exist is very short, about 10~^ s. A real observation sample 



^This corresponds to the time scale of the inner gap sparking 
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with extremely high time resolution (e.g. a few //s) will detect one point from only one transient 
"subpulse". In this case, a strong (even near 100%) circular polarization could be observed. 
Observational evidence was presented by Cognard et al. (1996). However, if sampling time is 
longer (e.g. ms), what one receive is the incoherent summation of emission from hundreds of 
transient beams at a given longitude. This largely diminishes the circular polarization, but does 
not affect much on the linear one. 

3.3. Polarization features of integrated pulse profiles 

Integrated profiles of pulsars are characterized by linear polarization, circular polarization, 
the sweep of position angle, and substantial amount of unpolarized emission. Our simulation 
shows these characters can all be reproduced. In our model, at every longitude, an observer can 
see the scattered emission from particles in the many sparking-produccd bunches around his line 
of sight. To get the integrated profile, we sum up incoherently the Stokes parameters of the 
scattered waves from all these bunches. For the simplicity of simulation, we assumed the particles 
in one bunch have the same Lorentz factor, since changing the Lorentz factor or considering the 
energy distribution of these particles should not alter the results following. 

We used the geometry presented in Fig. 5 for our model calculation, which helped to determine 
the different values for a sparking to every bunch in the ICS process. The key factor is the 

distribution of sparking spots in the polar cap. In the inner gap, a sparking is triggered when a 
7-photon encounters considerable B±, namely the perpendicular component of magnetic field, and 
produces an pair. A pair production cascade forms and the inner gap is discharged. This seems 
more likely to occur in a peripheral parts of the polar cap due to large B± while the central parts 
are inconducive to discharge. But at the edge of the gap, there is null electric field and thus it is 
impossible for a sparking to occur (Ruderman &; Suthland 1975). Therefore, the preferred region 
of gap discharge is located between the center and the edge on the polar cap. We assume in our 
simulation that the sparking spots have a gaussian distribution radially from the magnetic axis, 
with a maximum at ^ ~ 0.75^pc, where is the radius of the polar cap. However, changing the 
form of this radial distribution will not affect our results qualitatively. 

In the inner gap model, the electric potential across the gap is produced via monopolar 
generation. When the magnetic axis inclines from the rotation axis, the potential which accelerates 
the particles are different in the part of gap nearer to the rotation axis compared to that in the 
further part. Therefore, the probability of polar gap discharging varies with the distance to the 
rotation axis. This distribution asymmetry with respect to the magnetic axis has been included in 
our simulations. The probability of the sparkings was assummed to decrease exponentially with 
the azimuthal angle from the projection of the rotational axis. It is several times larger in the 
nearest part to the rotation axis than that in the farthest part. 

In Fig.8, we present a set of polarization profiles with various impact angle /3. We found that 
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the integrated profiles are highly hnear polarized in general, and their polarization angles follow 
nice S-shapes. One can immediately see that the position angle curves have a larger "maximum 
sweep rate" for a smaller /?, in excellent accordance with the rotating vector model. 

One important result is the antisymmetric circular polarization found in the central or core 
component, which can be up to 20%. As we mentioned before, circular polarization can be up 
to 100% on some part of "subpulses", but depolarization occurs in the process of incoherent 
summation of circular polarization of opposite senses. This also produces a substantial amount of 
unpolarized emission. For the linearly polarized emission, position angles of transient "subpulses" 
vary, causing further depolarization. 

However, in the conal components, circular polarization is insignificant, mainly because it is 
originally weak in the minibeam (see Sect. 3. 2). Since there is negligible variation of position angles 
of "subpulse" , almost no linear depolarization presents, and therefore, the conal components are 
always highly linearly polarized. 

4. Discussion and Conclusions 

Qiao (1988) and Qiao & Lin (1998) proposed an inverse Compton scattering model for radio 
pulsar emission. The simulations presented in this paper show that linear as well as circular 
polarization can be obtained from coherent inverse Compton scattering process from a bunch 
of particles, although the radio radiation of a single electron is highly linearly polarized. In the 
simulation process, there are two assumptions we have made, namely the coherence of the outgoing 
waves and small velocity dispersion of particles in a bunch. The circular polarization is favorable 
to be produced in emission regions near the polar cap in the ICS model. Observational data show 
that the 'core' emission is radiated at a place relatively close to the surface of the neutron star 
(Rankin 1993), the 'inner cone' is emitted in a higher region, and the 'outer cone' in the highest. 
Our simulations show that circular polarization tends to appear in the core components, while the 
'outer cone' emission is linearly polarized. 

One might have noticed that the circular polarization from our simulation of the ICS model is 
always strong and antisymmetric in the central part of integrated pulsar profiles. It is of intrinsic 

origin. However, almost no circular polarization appears in the conal components according to the 
ICS model, therefore the observed one must be converted from the strong linear polarization in 
propagation process (Radhakrishnan &; Rankin 1990; Han et al. 1998; von Hoensbroech & Lesch 
1999). 

One important ingredient to produce the circular polarization in the ICS model is the 
non-radial asymmetry of sparking distribution. We found from our summation that, without 
the non-radial asymmetry, circular polarization will be canceled by incoherent summation. The 
same conclusion has been reached by Radhakrishnan &; Rankin (1990) when they studied circular 
polarization using curvature radiation. 
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The polarization of the individual pulses or micro-pulses, however, should be related to the 
time resolution of observations according to the ICS model. With high time resolution (e.g. short 
to ~ lO^s), each sample will take emission from few transient beam that has strong linear and 
circular polarization. Depolarization happens during the sample time of a low resolution. This 
can serve as a test for our model. 

We emphasize that the width of observed subpulses is not necessary associated with I/7 
suggested by Radhakrishnan Sz Rankin (1990). The width of real subpulses observed from 
pulsars is generally 2-3 degrees, which was used to estimate 7 to be less than several tens. In 
our simulation, the width of transient "subpulse" produced by a bunch of particles along certain 
field lines reflect the width of the minibeam, which is mainly determined by the dimension of the 
bunch. The bunch radius near the star surface was taken to be 10 meters, corresponding to about 
one degree in longitude for core emission, and a little bit larger for cones. The width of a real 
observed sub-pulse might be related to the dimension of magnetic tube for the bunches more often 
to flow-out. Therefore, the upper limit of 7 can be released to 10^ or even higher. 

In our simulation, the low-frequency waves are assumed to be monochromatic. We have also 
made simulations for various a;in, and found that the polarization properties of transient beam 
and integrated profiles are quite similar to the results we presented above. 

We are very grateful to Dr. Bing Zhang and B.H. Hong for many extensive discussions, to 
Drs. R.N. Manchester, Q.H. Luo, Bing Zhang and anonymous referee for very helpful comments 
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Appendix: The amplitude of scattered wave in strong magnetic field 

We consider the complex amplitude of scattered wave in the electron-rest-frame, and then 
make Lorentz transformation for our case. The equation of motion for an electron with 7 = 1 in a 
magnetic field B is 

m r= eB{t) + xB. (Al) 

Here, e is the electron charge, and c is the speed of light. For conditions |r| -C c and Win ^ ^c, the 
magnetic component of the incident wave and the radiation damping force have been neglected. 
The electric field of the incident wave is taken as E(t) = Eq exp(— it^int) = e"^ E'q exp(— iwint) 
(here e™ is the polarization vector for incident wave) and the position vector of the electron as 
r{t) = roexp(— iwini). Generally, Eq and ro are complex vectors. One then can get 

ma;f„ro + eEo - i^ro x B = 0. {A.2) 
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For the incident wave in x-z plane (see Fig.l), fundamental linear polarization vectors 
of incident wave were defined e™ = {— cos^in, 0, sin^in} and — {0; 0}- The 

magnetic field B = {0, 0, B}. The fundamental vectors of scattered waves were chosen to be 
^i"* = {— cos 6*001 cos ^, — cos 6*001 sin ^, sin^out}) and 62"* = {sin^, — cos0, 0}. See Fig. 1 for the 
parameters of geometry. For a given incident wave, e™ = cosr/e™ + sin 7762'. Here e"^ is a complex 
vector, and hence 77 is a complex number generally. One can get the solutin of Eq.(Al) as: 



ro = — '—^ I sin — i/Li cos cos | , (^-3) 
^ cos?7sin0in(/x^ — 1) 



cos Tj cos 9i^ + ifi sin r] 



^2 



here jj. = Uc/i^m, v = Eqc^ / {eujf^) , and = In case of strong magnetic field, /x S> 1, and if 
■f] 7^ 7r/2, Eq.(A.3) can be simplified as being 

J'o = {0, 0, — i^rc COST/ sin 0}. (^-4) 

The electric field for the scattered wave can be written as (for |r| -C c) 



Es(0 = ^[n X (nx r)]i.ct 

= cos t] sin 6*111 sin 6'out exp[i^ (n • D - ct)\ ef"*. 



(^.5) 



Here, n = {sin ^out cos sin ^out sin (/>, cos^out} and D is the position vector from a scattering 
point to an observer, D = |D[. From Eq.(A.5) one sec that the scattered wave has only the term 
of 6°"*, and there is no 62"* term, which implies the completely linear polarization of the scattered 
waves. The differential cross section is 

here dQout is the solid angle of the scattered wave. The total cross section is 

= J dcr^ = (TThCOS^ 77 sin^^in, (^-7) 

here, fiTh = ^'''c is the Thomson section. For rj = (i.e. just for one polarization of e™), one 
can get cr|jy=o = 'T'tii sin^ ^in, which is consistent with the Eq.(16) at the lower frequency limit 
(wc > t^in) in Hcrold's (1979) when ^^j^ = 0. 

In the extreme case oi 7] = tt/2, the value of |ro| is extremely small (with a factor of (win/i^c)^) 
for the scattered waves in the radio band, so we will not discuss it in following. 

For an electron with 7 > 1, one can easily found through Lorentz transformation that the 
polarization vector of a wave is the same in an alternative inertial frame, but the complex amplitude 
of electric field of the wave in moving frame is 7(1 — /?cos^in) times that in a rest frame. In the 
electron rest frame (the frame moving with a electron with Lorentz factor 7), the complex amplitude 
of the electric field for the incident wave is Eq = 7(1 — (3cos9in)Eo. The complex amplitude for 
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the scattered wave in the electron rest frame is = — ^£^0 cos rysin^in sin0out/[7(l — /^cos^out)]) 
hence, in the laboratory frame it becomes 

E, = 7(1 + /? cos Ct)^f = -^Eo cos 7] sin sin Oout [7(1 - /3 cos ^out)]"^- (^-8) 

Therefore, for a moving electron with Lorentz factor 7, the electric field of the scattered wave 
becomes 

E,(t) = -^£;ocosr7sin^insin0out [7(1 - /^cos^out)]"^ exp[i— (n • D - ct)] e^*- (^-9) 
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Fig. 1. — The geometry of the ICS process for a single electron. The incident photon is coming as 
cjin in the x-z plane, and the scattered photons are going out around the magnetic field line (i.e., 
the z axis). The geometrical parameters used in the rest of the paper have been marked. In x-z 
plane, = 0. 



12 3 

Ooui (in 1/7) 

Fig. 2. — The total and linearly polarized intensities vary with the distance from the field direction, 
forming an micro-cone of out-going photons. The curves were calculated for = 30° and 80°. 
Other parameters were taken as: B = O.lBq (here Uc = 7.75 10^^ 7 = 100, Uin = 10^ s~^. 




Fig. 3. — The resonantly scattered waves: the total, linear and circular polarization intensity 
varying with the distance from the magnetic field direction were plotted as a solid line, dashed line 
and dash-dot-dash line, respectively. 
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Fig. 4. — The frequency dependence of total intensity (lower plots) and circular polarization 
percentage (upper plots) of the scattered waves about the resonant frequency, in four given directions 
of the out-going photon beam as (0,6lout) = (0°,0.5/7), (45°,0.5/7), (0°,2.0/7), and (45°,2.0/7). 
See Fig. 1 for geometry. The solid lines are plotted for the azimuthal angle ^ = 0°, and dash-dot- 
dot-dot lines for (f) = 45°. The left is plotted for emission at the angular radii ^out = 0-5/7, ^^'^ the 
right for 6lout = 2/7. 
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Fig. 5. — The geometry of the coherent ICS process. The sparking occurs at a point S above the 
polar cap, and produces the incident waves that are scattered by a bunch of particles. The emission 
visible to an observer at no is from a ring-like area satisfying • ng ~ cos(0.6/7). Point A is on 
this ring. 
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Fig. 6. — Polarization features of transient (a) core and (b) inner conal minibeams. The vectors 
stand for the polarization vector, their length is proportional to intensity. The filled and open circles 
represent the circular polarization of two senses, their size is proportional to circular polarization 
percentage. The circular polarization at two points in the core minibeam marked by the crosses 
gets to the maximum of 100%, with position angle jumps and null linear polarization. The sizes of 
circle in Fig. 6b have been amplified by 4 times for clarity. The thick line is the projection of the 
curved magnetic field. The line of sight goes across the beam, for example, along lines a, b and c 
in Fig. 6a and d in Fig. 6b. The radius of bunch is i?b = 10 m, magnetic inclination angle a = 45°, 
the bunch location (j)c = 225°, and 70 = 2000, cJout = 1 GHz. 
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Fig. 7. — Cross-sections of transient beams, cut by the lines of sight a, b, c and d in Fig 6. The 
obvious features of these "subpulses" are the presence of different circular polarization patterns and 
the variation of position angles. The circular polarization percentage and position angle variation 
range of inner cone emission (d) both are much less than that of core emission (a, b, c). 
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Fig. 8. — A set of integrated pulse profiles, simulated for the normalized impact angle (5^ = 0.15, 
0.30, 0.45, 0.60, 0.75, and 0.9. An S-shaped position angle swing, linear polarization, an 
antisymmetric type of net circular polarization and substantial amount of unpolarized emission 
can be presented from the ICS model. A small beam-shift between core and conal components due 
to different emission heights has been included in the simulation. 



